Statistically interacting quasiparticles in Ising chains 



Ping Lu 1 , Jared Vanasse 1 , Christopher Piecuch 1 , Michael 
Karbach 1,2 , and Gerhard Miiller 1 

1 Dcpartment of Physics, University of Rhode Island, Kingston RI 02881, USA 
2 Fachbereich Physik, Bergische Universitat Wuppertal, 42097 Wuppertal, 
Germany 

PACS numbers: 75.10.-b 

Abstract. The exclusion statistics of two complementary sets of quasiparticles, 
generated from opposite ends of the spectrum, are identified for Ising chains with 
spin s = 1/2,1. In the s = 1/2 case the two sets are antiferromagnetic domain 
walls (solitons) and ferromagnetic domains (strings). In the s = 1 case they are 
soliton pairs and nested strings, respectively. The Ising model is equivalent to a 
system of two species of solitons for s = 1/2 and to a system of six species of soliton 
pairs for s = 1. Solitons exist on single bonds but soliton pairs may be spread 
across many bonds. The thermodynamics of a system of domains spanning up to 
M lattice sites is amenable to exact analysis and shown to become equivalent, in 
the limit M — > oo, to the thermodynamics of the s = 1/2 Ising chain. A relation 
is presented between the solitons in the Ising limit and the spinons in the XX 
limit of the s = 1/2 XX Z chain. 



1. Introduction 

Ising chains are among the simplest systems of interacting degrees of freedom and have 
been thoroughly studied in a wide variety of circumstances including the presence of 
transverse fields, time-dependent fields, inhomogeneities in field or coupling etc. Is 
there anything of substance left that we can still learn from the Ising model in one 
dimension lj with homogeneous coupling? This paper presents a case (by no means 
the only one [2 [3]) for an affirmative answer. 

The Hamiltonian of the spin-s Ising model for s = 1/2, 1,3/2, ... on a periodic 
chain of N sites reads 



where the exchange coupling is antiferromagnetic (ferromagnetic) for J > ( J < 0) 
and h is a magnetic field. This model system has simple product eigenstates, a 
dispersionless spectrum, and no intrinsic dynamics to speak of. Its thermodynamics, 
derived via transfer matrix [H[S], is predictably simple. 

One interesting aspect of Tl s that promises usefulness in a wider context is the 
quasiparticle composition of the product eigenstates as will be demonstrated. The 
entire spectrum of Tl s can be systematically generated from opposite ends by different 
sets of quasiparticles with exotic exclusion statistics. The nature of these quasiparticles 
strongly varies with s but a systematics in their make-up is recognizable. 



N 




(1.1) 
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In H1/2 we consider antiferromagnetic domain walls (solitons with spin ±1/2 and 
fractional exclusion statistics) for J > or with ferromagnetic domains (strings of 
flipped spins with integer- valued exclusion statistics) for J < 0. The corresponding 
quasiparticles in Hi turn out to be soliton pairs with spin 0, ±1 (for J > 0) and nested 
strings (for J < 0), both with unusual exclusion statistics. 

We use the concept of statistically interacting quasiparticles to show that the 
thermodynamics of H s is equivalent to that of a gas of solitons (for s = 1/2) or soliton 
pairs (for s = 1). The same framework is shown to work also for the thermodynamics of 
string particles. It is expected that these particles, whose detailed exclusion statistics 
is worked out here, are still relevant in integrable spin chain models away from their 
Ising limit. The particles identified here then become objects of a coordinate Bcthc 
ansatz [6J El [8] applied to those models. 

We first review the concept of statistical interaction and its use in thermodynamics 
(Sec. HJ. Then we introduce the soliton particles for H1/2, describe their exclusion 
statistics, and determine the thermodynamics in a magnetic field from a soliton 
perspective (Sec. [3]). Next we introduce the six species of soliton-pair particles that 
govern the spectrum of Hi and work out their thermodynamics in zero magnetic 
field (Sec.[4]). Then we present the combinatorics for the statistical interaction of a 
system of strings in H1/2 and of nested strings in Hi. We proceed by calculating 
the thermodynamics of a system of strings of restricted size and recover the Ising 
thermodynamics when that restriction is lifted (Sec. [5]). Finally, we assess the progress 
reported here and discuss possible extensions and comparisons (Sec. |6|) including a 
relation between solitons and spinons, both with semionic statistics ( | Appendix A[ ). 

2. Statistical interaction 

Quasiparticles in solid matter are not restricted to be either bosons or fermions. In 
integrable quantum many-body model systems (6j El H] quasiparticles with infinite 
lifetimes and unusual exclusion statistics have indeed been identified. The generalized 
Pauli principle as introduced by Haldane [9] expresses how the number of states 
available to one particle is affected by the presence of other particles: 

A,/, • ^ry,,A.V,. (2.1) 

3 

The indices i,j refer to distinct particle species. The gij are statistical interaction 
coefficients. For bosons we have = and for fermions g^ — Sij. Upon integration 
Eq. (|2.ip becomes 

di=Ai-Y,m{Nj-5ij), (2.2) 
3 

where the A4 are statistical capacity constants. The number of many-body states 
containing {Ni} particles of the various species is then determined by the multiplicity 
expression 

wmn - n ( dl + n* " 1 ) ■ ( 2 - 3 ) 

We shall determine the ingredients Ai,gij to (|2.2[) for two species of solitons or 
N species of strings, all pertaining to H1/2, and for six species of soliton pairs or 
N(N + l)/2 species of nested strings pertaining to Hi. 
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The thermodynamic properties of a macroscopic system of statistically interacting 
particles are amenable to a rigorous analysis as shown by Wu [TO]. For given sets of 
1-particle energies e,;, statistical interaction coefficients gij, and statistical capacity 
constants A4, the grand partition function is 

i 

where the quantities Wi are determined by the nonlinear algebraic equations 

^ = Mi+^-I>a»(i±f). (2.5) 

The temperature T and the chemical potential fj, are the control variables. Additional 
control variables such as external fields may come into play as part of the energies Cj. 
The average numbers of particles, (Ni), of each species are related to the Wi by the 
linear equations 

Wi{Ni)+^gij{Nj) = Ai. (2.6) 
3 

We shall apply this method of exact analysis to the solitons and the strings in the 
context of TL1/2 and to soliton pairs in the context of TL\. 

3. Solitons 

Here we consider TL1/2 with J > and h > for even or odd N. The task at hand 
has a combinatorial part and a statistical mechanical part. We first relate the Ising 
spectrum to soliton configurations, then we undertake a thermodynamic analysis of 
the soliton system, using the methodology outlined in Sec. [2] 

3.1. Combinatorics of solitons 

Among the four distinct bonds in the general product state \<j\<Ji ■ ■ ■ <t/v) ( see Tablc[T]), 
the bonds ff, JJ represent solitons with spin +1/2, —1/2, respectively, and "fj,, J,| arc 
vacuum bonds. Close-packed solitons with like spin orientation reside on successive 
bonds (e.g. TTT)i whereas close-packed solitons with opposite spin orientation are 
separated by one vacuum bond (e.g. TTTT)- More generally, the number of vacuum 
bonds between nearest-neighbor solitons with like (opposite) spin orientation is even 
(odd) . Solitons only interact statistically. The energy of a soliton is unaffected by the 
presence of other solitons. 



1 + Wi 



(2.4) 



Table 1. Distinct bonds in H1/2, their soliton content, and their contribution to 
the energy of the product eigenstate (relative to the soliton vacuum). 



bond 


TT 


ii 


TT 


TT 


N+ 


1 
















1 








AE 


J+h 
2 


J-h 
2 
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The two soliton vacuum states, | T-L • " TJ.) an d | IT • ■ ■ IT), represent the lowest 
energy level for even N. The lowest level for odd N is 27V-fold degenerate and contains 
one soliton. The soliton content of an Ising eigenstate is specified either by the numbers 
N± of spin- up/down solitons or, alternatively, by the total number of solitons and the 
magnetisation: 

N A = N + + N_, M Z = ^(N + -N_). (3.1) 

The energy level of all states with Na solitons and magnetisation M z is 

E{N A ,M Z ) -E = ^N A J + hM z , (3.2) 

where Eq = —NJ/4 is the energy of the soliton vacuum. 

How many Ising eigenstates exist for given N + and iV_ (or Na and M z )l 
The solution of this combinatorial problem is the following multiplicity expression 
constructed from extensive tabulated data such as sampled in Table [2] 

^♦."-'-j^nC^*- 1 ). m 

a— i ' 

d„ = \{N-l)- l -Y^ N °'-^-\ (3-4) 

a' 

It is indeed compatible with the standard form (|2.3[) . The range of Na is 0, 2, • • • , N 
for even N and 1,3, • • • ,N for odd N. Na = N is only realized for the two states 
with N + = Na or iV_ = Na- This multiplicity expression specifies the statistical 
interaction between soliton particles. 



Table 2. Number of states, W A (N + ,N_), with N A = N+ + N_ solitons and 
magnetisation M z = (N+ - N-) /2 for H 1/2 with N = 6 (left) and N = 7 (right). 



M Z \N A 





2 


4 


6 




3 








1 


1 


2 






G 




6 


1 




9 


G 




15 





2 


12 


G 




20 


-1 




9 


6 




15 


-2 






G 




6 


-3 








1 


1 




2 


30 


30 


2 


64 



M Z \N A 


1 


3 


5 


7 




7/2 








1 


1 


5/2 






7 




7 


3/2 




14 


7 




21 


1/2 


7 


21 


7 




35 


-1/2 


7 


21 


7 




35 


-3/2 




14 


7 




21 


-5/2 






7 




7 


-7/2 








1 


1 




14 


70 


42 


2 


128 



3.2. Statistical mechanics of solitons 

For the statistical mechanical analysis of H1/2 as a soliton gas we use the statistical 
capacity constants A a = (N — l)/2 and the statistical interaction coefficients g aa i = 
1/2 extracted from Eq. (13. 4j) . and the soliton energies e CT = ( J + <rh)/2 from TablcQ] 
We have to solve two coupled nonlinear algebraic equations of the type (12. 5|) : 

=nl + w± ) + -ln — + -ln ^— . (3.5) 

2fc B T 1 ±7 2 l + w± 2 l + w T V ; 
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The solution, 



W± = 7: 



determines the grand partition function via 
for the capacity constants. The result, 



±h/k B T _ 1 + ^J( e ±h/k B T _ iy + 4:( ,(j±h)/k B T 

with the asymptotic value A± 



J1K 



coshff 



V sinh' 2 



H 



N 



K 



J 



H 



(3.6) 

> N/2 

(3.7) 



4fc B T' 2fc B T' 
is in exact agreement with the well-known canonical partition function Zn obtained 
via transfer-matrix 0]. In the relation Z = Zpje NK , the factor e NK accounts for the 
relative shift of energy scales used in the two methods. 

For the average numbers of solitons, (N±), we infer from (|2.6p the two coupled 
linear equations, 



w± 



1, v N 



which have the solutions 



±H 



Vsinh 2 H + e~ iK ± sinh H 



N/2 



-IK 



(3.8) 



(3.9) 



sinh 2 H + e~ 4K + cosh i/Vsinh 2 H 
The dependence of (N + )/N on J/fc^T is shown in Fig.[T]for various values of h/J. All 
curves start from (N+) /N — 1/4 in the high-T limit. For ft = Owe have (N+) = (N-); 
this curve has a monotonically decreasing trend toward zero as T — > 0. For h > 
(ft < 0) the average number (N + ) of solitons with spin directed antiparallel (parallel) 
to ft, is more (less) rapidly suppressed as T — > 0. For sufficiently weak, negative fields, 
> h/J > —0.25, the curve is still monotonically decreasing. For —0.25 > ft/J > — 1, 
it acquires a smooth maximum at finite, nonzero T. For ft/J < —1 the curve is 
monotonically increasing toward (N + ) /N=l. Here the ground state contains N spin- 
polarized solitons. 




Figure 1. Average number (N+)/N of spin-up solitons per lattice bond versus 
J/ksT for several values of magnetic field. 
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4. Soliton pairs 

Here we consider Tli with J > and /i = for even or odd N. The generalization 
to h 7^ is straightforward conceptually. The independent particles are now soliton 
pairs. 

4-1- Combinatorics of soliton pairs 

The nine different kinds of bonds are listed in Table [31 Each bond can accommodate 
up to two solitons. The energy of a soliton is not the same in all configurations. 

Table 3. Distinct bonds in Hi, their soliton content, and their contribution to 
the energy of a product eigenstate (relative to the soliton vacuum). 



bond 


TT 


oo 


41 


T° 


°T 


1° 


°i 


Tl 


IT 


N+ 


2 


1 





l 


l 














iV_ 





1 


2 








1 


1 








AE/J 


2 


1 


2 


l 


l 


1 


1 









Our search for the independent particles again starts from extensive tabulated 
data for Wa(N + ,N-) such as sampled in Table H] Several clues suggest that the 
independent particles are soliton pairs. For example, the number of solitons is always 
even. Also, the number of states with N- = grows cx N N +/ 2 for N + <C N as opposed 
to the growth cx N N + observed in H\/ 2 

The systematic examination of the data tables for Wa(N+, N-) points to the 
existence of six distinct species of soliton-pair particles, two groups of three species 
with spin (+1,0, —1). In the first group the paired solitons are confined to the same 
bond. In the second group the paired solitons are deconfined. They can be on bonds 
with any number of lattice sites between them. 

Confined-soliton pairs with spin up (named r+) are identified by any element | 
in the product state. In like manner, confined-soliton pairs with spin zero (down) are 



Table 4. Number of states, Wa(N+, N-), with N A = N+ + N- solitons and 
magnetisation M z = (N+ - N_) /2 for Hi with N = 3 (left) and N = 4 (right). 







4 






M Z \N A 





2 


4 


6 


8 




M z \Na 


2 


6 




4 










1 


1 


3 






1 


1 


3 








4 




4 


2 




3 




3 


2 






6 


4 




10 


1 


3 


3 




6 


1 




4 


8 


4 




16 





6 




1 


7 





2 




16 




1 


19 


-1 


3 


3 




6 


-1 




4 


8 


4 




16 


-2 




3 




3 


-2 






6 


4 




10 


-3 






1 


1 


-3 








4 




4 




12 


12 


3 


27 


-4 










1 


1 












2 


8 


48 


20 


3 


81 
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named r (r_) and identified by elements oo (J.J,) in the product state. Deconfined- 
soliton pairs with spin up (down) are named q + and identified by any element 
f o ■ • • o ] ([ o ■ • ■ o I) in the product state, where the presence of n — 1,2,... site 
variables o between two site variables T (J.) indicate the presence of n — 1 spin-zero 
confined-soliton pairs (?"o). Deconfined-soliton pairs with spin zero are named qo and 
are identified by elements |o---o|or|o---o|in the product state. 

A list of name, motif, and soliton content for all six species of soliton-pair particles 
is shown in the top three rows of TableO In some instances, two close-packed particles 
share one lattice site (e.g. o o o, f o ff, | o J, o J,), in other instances, there is one 
vacuum bond in between (e.g. TTI-l> T ° Tl ° ?)• The particle ro can only exist inside 
one of the particles q+, qo, g_. The former is instrumental to the soliton deconfinement 
in the latter. 

The six particles are thus naturally classified into three groups, the confined- 
soliton pairs r+, r_, the deconfined-soliton pairs q+,qo,q-, and the spacer particle ro 
(deconfinement agent). In a vague QCD analogy, solitons play the role of quarks, the 
soliton pairs r+, r_, g+, <7oj 9- the role of mesons and baryons, and the spacer particle 
tq the role of gluon with opposite action. 



Table 5. Specifications of particles in "H\: confined-soliton pairs (r_|_,r_), spacer 
particle (ro), deconfined-soliton pairs (q+, qo, <?— ); motif in product state; soliton 
content; index m used in 14.lt ; statistical capacity constants A m ; energies e m . 



particle 


r + r_ 


7"0 


9+ 9- qo 


motif 

N+ + iV_ 


TT U 

2+0 0+2 


OO 

1 + 1 


T ° T l°l T°l, l°T 

2 + + 2 1 + 1 


m 


1 2 

JV-l AT-l 

2 2 

2J 2J 


3 

J 


4 5 6 
f-1 f-1 AT-2 
2J 2J 2J 



We have determined the multiplicity expression 

^j)=^n(*" + i™- 1 ). («> 

m— 1 x 7 

^?7i ^-m ^ ^ 9 mm' {^Cm 1 ~~ ^mm')' (^ - ^) 
rn ' 

iV E = Xi + X 2 + Z 3 + 2(X 4 + X 5 + X 6 ) < AT, (4.3) 

for product eigenstates containing X m soliton pairs of species m — 1, . . . , 6, where the 
index m is defined in Table [5J It confirms the independent status of the soliton-pair 
particles and contains the specifications of their statistical interaction. The capacity 
constants A m and the particle energies e m are given in Table [51 and the interaction 
coefficients g mm ' in Table [5] Again there exist restrictions and exceptions regarding 
the allowed configurations {X m }. We do not list them here because they have no 
bearing on the statistical mechanical analysis. The only model specifications needed 
are the quantities A m , e m , g m m'- 
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Table 6. Statistical interaction coefficients g mm i between soliton-pair particles 
as identified in Table [5] 



Qmm' 


1 


2 


3 


4 


5 


6 


l 


l 

2 


l 

2 


i 

2 





1 


i 

2 


2 


1 


1 


1 


1 





1 




2 


2 


2 


2 


3 











-1 


-1 


-1 


4 


1 


1 


1 


1 


1 


1 




2 


2 


2 


2 


5 


1 

2 


1 

2 


1 

2 


1 


1 


1 

2 


6 


1 


1 


1 


2 


2 


2 



Note that the g mm > include some zeros and some negative values. To make sense 
of these peculiarities we rewrite each of the six binomial factors of (|4.ip in the form 



(4.4) 



where 

— A m -f- gmmi Y m — ^ ^ gmm' X m ' ~t~ 1- 0^'^) 

m'^m 

The maximum capacity for particles of species m, 

X m = — . (4.6) 

9mm 

is thus primarily dictated by the diagonal coefficient g mm , but is also influenced by 
the off-diagonal coefficients g m m' via Y m . 

If one of the off-diagonal coefficients is zero, gmm' — for m! ^ m, it merely 
means that the presence of particles of species to' has no effect on the capacity for 
particles of species m. If one of the diagonal coefficients vanishes, g mm = 0, then (|4.4[) 
does no longer limit the capacity for particles of species m. This can either mean that 
there is no limit (as is the case for bosons) or it can mean (as is the case here for 
m = 3) that a limit is implied by a different rule associated with (|4.1|) . 

The existence of negative off-diagonal coefficients g mm ' a s found in Table [S] for 
m = 3 and m! = 4, 5, 6 has the consequence that particles from species m! contribute 
negatively to Y m . Adding particles of species m! increases the capacity of the system 
for particles of species to. This is indeed to be expected because the latter can only 
exist inside the former. In this instance, the Pauli exclusion principle turns into what 
might be called an accommodation principle. 

4-2. Statistical mechanics of soliton pairs 

Carrying out the statistical mechanical analysis of 7ii as a gas of soliton pairs starts 
with solving six coupled nonlinear equations of the type (|2.5|) with /i = 0, the e rn from 
Table El and the g mm ' from Table [H Symmetry implies wi — W2, — w^. The 
remaining four equations in exponentiated form (with K = —J/ksT), 

_ 2K _ WiW4W 6 _ K _ (l+W 3 )wiW4W 6 

~ (1 + W4)(1 + W 6 )' 6 ~ (1 + Wi)(l + W4)(l + W 6 )' l ' j 
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-2K _ Wl(l + W 3 )w|w| _ 2K _ Wi(l + W 3 )w 4 wl 

~ (1 + W!)W 3 (1 + W 4 )(l + W & ) 21 6 ~ (1 + Wl)w 3 (l + W 4 )(l + W 6 ) : 



can be simplified into 

1 + W3 K w 3 K ■, , 2 + w 6 2K 

=e , — =e , W4 = l + ui6, =e , (4.9) 

1 + wi We w±we 

and reduced to a quadratic equation for w 3 with the (physically significant) solution 




w 3 = cosh K- - + \ cosh if-- +2. (4.10) 



The grand partition function (I2.4|) with the (asymptotic) A m from Table [5] and the 
solutions (|4.9p . (|4.10p becomes 

wi) 2 (1 + w 4 ) 2 (1 + we) 2 1 ^ _ r n s jq Af _ 7 ™ ain 
2 Z2 -[(l+w 3 )e J - Z N e , (4.11) 

1 w i w 6 

in agreement with the transfer-matrix result for the canonical partition function Zn 
[5], where the factor e KN again originates from a shift in energy scale. 

For the average numbers of soliton pairs, (N m ), we must solve six linear equations 
of the type (|23|) with the w m from (jH)]) . |[4TTD]) . Symmetry dictates that (iVi) = (AT 2 ) 
and (iV 4 ) = (N 5 ). The solution reads 

iV w 3 {wl+2e K ) 
2 (w 3 + ^(wfe-^ + Aw 3 + 2 - 2e K ) : 



(N \- N 2(w 3 + l~e K ) 

[ 3> (w 3 + l)(w 2 3 e- K +4w 3 +2-2e K Y [ ' 

(N 4 ) = (N 5 ) = ^w 3 (N 3 ), (N 6 ) = ~w 3 (N 3 ). (4.14) 

The relation (Nq) = 2(N 4 ) = 2(N§) may be anticipated on the basis of the motif 
shown in Table[5] The reduced averages (N m )/N are plotted versus J/fc^T in Fig. [21 

Increasing the temperature from T = results in a gradual increase of average 
particle numbers from all species. Notice that the spacer particles tq (m — 3), which 
can only exist inside particles of species q+,qo, g_ (m = 4, 5, 6) are slowest to appear 
in significant numbers as would be expected. 

For the generalization of these results to h ^ we must add the Zeeman 
contribution to the energies e m listed in Tabled The statistical mechanical analysis 
of soliton pairs as demonstrated here is by no means limited to the Ising Hamiltonian 
(jl.ip . We can freeze out some of the particle species by making their activation energies 
infinitely large, e m — > oo. This has the consequence that w m — > and (N m ) — > 0. For 
the remaining active particles we can assign arbitrary values e m for their energies. For 
example, if we freeze out all particles except those of species m = 1, 2 then the results 
of Sec. [3] for Hi/2 are, effectively, recovered. 

Among the issues that must be heeded in generalizations of the calculations 
reported here to models with arbitrary particle energies are the following: (i) The 
particle ro can only exist inside a particle qo, q + , or Therefore freezing out the 
latter three will freeze out the former even if its energy remains finite, (ii) The particle 
go comes in two parity-violating versions. In the context of Tii or other models where 
the two configurations t ° 4 an d 4 ° T have the same energy, they can be treated as 
identical particles. However, in situations where the two configurations have to be 
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Figure 2. Average numbers (N m )/N of soliton-pair particles (per site) versus 
inverse temperature J/kgT of Hi. 



assigned different energies we must treat them as belonging to different species and 
determine their statistical interaction with each other and with all the other particle 
species, (hi) A spin interaction beyond nearest neighbors added to 7ii will, in general, 
produce a coupling between the particles listed in Table 



5. Domains and nested domains 



Here we consider Ti s for s = 1/2, 1, J < 0, and h < 0. In the following we describe 
how the entire spectrum is systematically generated from the ferromagnetic ground 
state | It • • " T) by domains (s = 1/2) or nested domains (s = 1) of flipped spins. 
These domains are independent particles subject to a statistical interaction. The 
thermodynamic analysis of domains is then carried out for the s = 1/2 case at h = 0. 



5.1. Combinatorics of domains 

In the notation used here, {\<ti • • • &N)}m represents the set of m product vectors with 
r flipped spins that are generated from \a\ ■ ■ ■ er/v) via translations. The 2 N = 16 
states for N = 4 in this representation are 

{1mm?, {muni {immi 

{\nnnl {iuu»i, {iuu»* (5.1) 

The first among them is the (non-degenerate) ground state of H-i/2 with J < and 
h < 0. Domains are strings of fi consecutive down-spins. In (|5.1[) the states in the 
second set contain one 1-string, and the state in the fourth set two 1-strings. The 
states in the third, fifth, and sixth set contain one string with ji = 2, 3, 4, respectively. 
Each string of length fi contributes the amount J + fih to the energy of the state. An 
Ising chain of length N can thus accommodate strings with [i = 1, . . . , N — 1, which 
are treated here as distinct species of independent particles. The lone state containing 
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one string with fi — N is exceptional in several respects, ignorable in macroscopic 
systems. 

What is the number of product eigenstates that contain a configuration {N^} 
of strings? Since there must be at least one up-spin between successive strings, only 
those configurations can be realized which satisfy the constraint 

N~l 

]T( M +1)A^ <iV. (5.2) 

We found that the number of states with given string configuration is determined by 
the multiplicity expression 

n(* + # _1 )- (5 " 3) 

JV-l 

d IM = N- f i-J2 <W'(^V - <W)i ( 5 - 4 ) 

M' = l 



where 



9nfi' 



( i i 



5.J2. Combinatorics of nested domains 

The concept of nested quasiparticles in lattice models is well-known in the context 
of the Bethe ansatz as applied, for example, to the Hubbard model or to integrable 
spin-1 models [H [HJ [12l H31 [14j [15l [16]. The nested particles in Ising product states 
have a particularly simple structure. 

In the context of Hi the nesting involves two shells. The particles on the outer 
shell (/x-strings) are structurally identical to the strings of Wi/2- We start from the 
/it-string vacuum, {| || ■ ■ • t)}?j an d generate a total of 2 N product states composed 
of site variables f and o. On the inner shell we take any /i-string of the outer shell 
and use it as the vacuum for ^-strings. Hence a i/-string is a sequence of v successive 
|-sites embedded in a region of \i o-sites between consecutive f-sites. Naturally, we 
must have v < p. This prescription is illustrated in Table [7] for N — 4. The two-shell 
nesting of string particles leads to the multiplicity expression 

with d^ from (|5.4|) . g^i, r from (|5.5p . and 

d ( + l) =n-v-Y,9 V u>{N { f -8 VV ,), t»=Y,vNW>. (5.7) 

u' v 

As in previous applications, there are instances (ignorable for macroscopic systems) 
where expression (|5 .6[) is inapplicable. 
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Table 7. Nested-string interpretation of the product eigenstates for N = 4 in 
H\. The 2^ states on the left represent the outer shell of the nesting. Each 
/i-string is underlined and serves as the vacuum for i/-strings on the inner shell. 
Prom each of the 2^ effective H1/2 product states on the left are thus generated 
one or several Hi product states on the right for a total of 3^. 
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5.3. Statistical mechanics of domains 

Returning to TL1/2 with J < and setting h — 0, we now derive the exact 
thermodynamics of a system of strings via the method outlined in Sec. [5J It is evident 
from Wu's derivation [TU] of Eqs. (|2.5p that their applicability in the present context is 
limited to situations where the system has a large capacity for strings of all sizes that 
are permitted. To circumnavigate this restriction we introduce a limit on the length 
of allowed strings, fj, < M <C N. The thermodynamic limit of H1/2 requires that we 
set N — > 00 before setting M — ► 00. 

With the specifications regarding statistical interaction of strings from Sec. 15.11 
we write for the grand potential the expression 

M , , 

n M (K) = -k B Ty / A f ,\n(^—), A„ = N-n, (5.8) 
where the satisfy 

4K = lnK + l)-^ P ^ln^^ ) K = (5.9) 

n'=l I" 

The transformation of variable, £ M = ln(u> M + 1), turns Eqs. (|5.8p and (|5.9p into 

I J\ M 

n M {K) = L±Y,( N -ri ln (i - e ^") > ( 5 - 10 ) 

M M 

^ = iK-fi ln(l-e-V) - ^2 ln(l -e"^') . (5.11) 
Introducing the quantity 

= (5.121 



4K ^ 
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we rewrite (|5.1ip in the form 

^ = 4X[1 + ( M +1)$ M -$„_!]. (5.13) 
This sets the stage for determining $m via a recursive scheme: 

<f> 1 = -±ln(l- q i+™"), q = e^ K , 

$2 = $1 - ^ In (1 - exp (-4X(1 + 3* M ) - In (l - 9 1+2 *«))) , 

leading to 

$„=- — ln(l- g 1+2 *" 7 T ). (5.14) 

Setting fi = M we arrive at a polynomial equation for g* M : 

qq {M+2)* M + (1 _ g ) g 2*M _ 2g *« +1 = 0. (5.15) 

The solution of (|5. 15|) substituted into (|5. 10|) via (|5 . 14[) and (15.13P determines the 
grand potential of a system of strings with maximum length M in a chain of N sites 
with M <C N. Taking the limit N — > oo while keeping M finite we have 

1 171 M 

lom(K) =lim-n M (K) = Ll^^i-e-^) = -\J\* M . (5.16) 

M =i 

If we now take the limit M — ► oo, the first term in Eq. (|5 . 1 5[) vanishes, and the 
solution, 

q^={l + ^r\ (5.17) 
substituted into (|5.16|) . yields 

Uoo(K) = -^\n(l + e- 2K ), (5.18) 

which is indeed the exact result for H1/2 with h = 0, J < and the string vacuum at 
the origin of the energy scale. 

The statistical mechanics of a system of domains with maximum length M on 
a lattice of N sites may very well be of interest in a number of contexts outside 
magnetism. We have reduced the problem to solving a polynomial equation of degree 
M + 1. Consider the entropy per site of strings with fi < M on an infinite lattice, 
sm{K) = limiv-^oo Sm{K)/N, inferred from (|5.16p . Compact analytic solutions are 
readily calculated for M = 1 (one-strings only) and M = 00 (all strings allowed): 



k B \ 2 J l + 4e- 4K + Vl + 4e- 4 * 

^ = Mi + ^'> + ^. ("0) 

Entropy curves for several M are shown in Fig. [3] As we relax the restriction 
on the length of permissible domains, the entropy at any given nonzero temperature 
becomes larger. The relative contribution of longer domains is larger at low T than 
at high T. All domains have the same energy. With T increasing, the longer domains 
tend to be crowded out by the shorter ones. As the restriction on length is lifted 
altogether, the Ising result (M — 00) is approached from below. 

The same type of analysis is applicable to any model with spin- 1/2 Ising 
product eigenstates and with arbitrary energy values e m , /i = 1,2,..., M assigned 
to the domains. The left-hand side of Eq. (|5.9p must then be replaced by 4-K^, 
Kn = £fi/4kBT. For 7ix/2 at ft, 7^ we must use e M = J + fih. 



Statistically interacting quasiparticles in Ising chains 



14 




k B T/|J| 

Figure 3. Entropy per site for N — » oo versus reduced temperature of a system 
of domains with maximum length M. The case M = oo represents T~L\i<i. 



5.4- Distribution of domains 

What is the relative frequency of occurrence of domains of size fj, for given maximum 
size M at temperature T in an infinite chain? To answer this question we adapt Wu's 
linear equations (|2.6[) to the situation at hand: 

M \i 

w„(n ll )+ 2 53 vV) =1, n = l,...,M, (5.21) 

where n M = N^/N and where we have ignored a contribution of 0(n/N) to the 
right-hand side by effectively taking the limit N — > oo, while keeping M finite. The 
quantities = e^ fI — 1 are known from the solution of (|5 . 1 5|) via (|5.14[) and (|5. 13[) . 

Here we carry out the calculation for the case M — > oo. The solution (|5 . 1 T|) 
substituted into (|5.14[) yields 

„*„ _ 1 (, , v 7 ? 



which, upon substitution in (|5.13[) . produces the needed in (|5.21[) : 

^ = 1 + (1 + yffi" (M = ^ (5 

Now we rewrite Eqs. (|5.21[) in the form 

«' /1 <n p )+53(n»') = C, M=l ) 2,... (5.24) 
where the quantity 

oo 

c = i-J2"M ( 5 - 25 ) 

i/=i 
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can be treated as a constant to be determined self-consistently at the end. The solution 
of Eqs. (|5.24p . obtained by induction, is 

p / oo \ _1 n 

yi \ ^— * in. I J- X in. -I- I 

and, after normalization, 



From (|5.12[) and (|5.26p we infer that — g*^ a quantity evaluated in (|5.22[) . The 
assembly of the ingredients (|5.22p . (|5.23[) . and (|5.17p to expression (|5.27p produces the 
following explicit result for the distribution of lengths /i of string particles in Tii/ 2 at 
temperature T and zero magnetic field: 

fa e~ 2K 
^ )= (TT^ = (l + e-^ ' M=1 ' 2 '- (5 ' 28) 

This is a realization of Pascal's distribution, P(fi) = 7(1 — 7) A1 ~ 1 , if we set 7 = 
e~ K I (e + e ). This result was previously derived by Denisov and Hanggi |3J using 
a very different method in their study of finite Ising chains with open boundaries. 
This distribution indeed favors short strings in the crowded conditions at high T, in 
agreement with observations made in our discussion of the entropy curves (Fig. [3]). At 
low T the distribution is flat, consistent with the fact that all strings have the same 
energy. With some additional effort our solution can be generalized to finite M, and 
to models with arbitrary values for the energies e M of domains of size fi. 



6. Conclusion 



We have demonstrated that the conceptual framework of statistical interaction 
between quasiparticles in many-body systems [SJ [TU] leads to significant new insights 
into the statistical mechanics of Ising chains and related models with spin- 1/2 or spin- 
1 product eigenstates on a one-dimensional lattice. We have identified, in particular, 
the nature of complementary sets of independent particles on the basis of which the 
spectrum of Ising chains with s — 1/2 and s — 1 can be generated systematically from 
either the ferro- or antiferromagnetic ground state. 

The Neel state is the vacuum for solitonic particles. In the s = 1/2 case the 
solitons themselves are the independent particles. They are antiferromagnetic domain 
walls, confined to single bonds, with spin ±1/2 and semionic statistical interaction. 
In the s — 1 case the solitons are merely building blocks of particles. All independent 
particles are soliton pairs. The paired solitons may be on the same bond or on bonds 
any number of lattice units apart. We have carried out the exact statistical mechanical 
analysis of solitons (two species) for s = 1/2 and of soliton pairs (six species) for s = 1. 

The state with all spins up is the vacuum for string particles. In the s = 1/2 
case the independent particles are domains of overturned spins and in the s = 1 case 
they are nested domains, i.e. domains inside domains of halfway overturned spins. By 
working out their exact statistical interaction we have set the stage for the statistical 
mechanical analysis of domains and nested domains. We have carried out that analysis 
for the s = 1/2 case and established contact with previous work based on different 
methods [3]. 
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The work presented here opens up numerous opportunities for extensions and 
comparisons including the following, (i) The methodology developed in Sees. |3] 
and 2] for the identification and specification of independent solitonic particles looks 
promising for applications to Ising chains with s > 1 and to Ising ladders. Preliminary 
results for W3/2, for example, indicate that the independent solitonic particles contain 
at least two and no more than six solitons. This again includes particles confined 
to one bond and particles spread across many bonds with more than one species of 
spacer particles acting as deconfinement agents, (ii) A question of considerable interest 
is how the methodology developed here can be generalized to situations with Ising 
interactions beyond nearest neighbors, which, in general, leads to a coupling between 
solitonic particles and between string particles, (iii) There exist integrable spin chain 
models with a parametric Ising limit. Consequently, the solitonic particles analyzed 
here must exist in some variant form away from the Ising limit of those models. One 
such link (to the spinons of the s = 1/2 XX Z model) is outlined in |Appendix A| 
Corresponding links are bound to exist between the string and nested-string particles 
of Ising chains and the string solutions of the Bethe ansatz applied to integrable spin 
models with s = 1/2 [HI [HI [HI [20l EH [22] and s = 1 pll 03 [13 El El QI] near 
their Ising limits. 

Appendix A. Solitons versus spinons 

The ground state of the s — 1/2 XX Z model, 

N 

T~txXZ = ^2 [j±(SnSn+l + ^n^n+l) + ^X^n^n+l] l (A-l) 

n=l 

at Jj_, J z > for even N is non-degenerate except in the Ising limit (J± = 0). The 
finite-size gap is of 0(N~ 1 ) in the planar regime ( J± > J z ) and of 0(e~ N ) in the axial 
regime (Jj_ < J z ). The lowest energy level in both regimes has been identified as the 
(unique) vacuum of spinons 22J. The two lowest levels, again in both regimes, can be 
identified as the (twofold) vacuum of solitonsQ] 

Spinons and solitons have similar but not identical semionic exclusion statistics. 
The similarities and differences are encoded in the multiplicity expressions. Equations 
(|3.3[) - (|3.4p for solitons are to be compared with 

W 8 (N + ,N-)=U ( da+ Nr 1 ), (A.2) 
da = \(N + l)-^(N <r ,-6 <ra ,) (A.3) 

a' 

for spinons [S]. Away from the Ising limit, solitons (and spinons) are dispersive 
and scatter off each other elastically. Both kinds of particles are most conveniently 
identified by their momentum quantum numbers. Every XX Z eigenstate has a unique 
spinon composition and a unique soliton composition. The relation between the spinon 
composition and the soliton composition is most transparent in the XX limit ( J z =0). 

In Ref. [23] a motif was developed that relates the configuration of (free) Jordan- 
Wigner fermions with the configuration of spinons. This motif is reproduced in Fig. lAll 
for N = 4 (16 eigenstates) and amended to also show the soliton configuration. The 

I The names attributed to quasiparticles in quantum spin chains vary among authors. Our usage is 
common but not universal. 
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allowed fermion momenta (in units of tt/N) are 

f {l,3,...,2iV-l} for even N F 
1 fc \ {0, 2, . . . , 27V - 2} for odd N F 1 ' 

and the allowed spinon orbital momenta (in units of ir/N) are 

mi = — , — +%.. .,N- —, (A.5) 

where Np is the number of fermions and N$ = N + + N_ the number of spinons in 
any given XX eigenstate. 

The exact spinon configuration is encoded in the fermion configuration as 
described in the following: (i) Consider the the gray fork as dividing the fermion 
momentum space into two domains, the inside and the outside. The outside domain 
wraps around at the extremes (m, = iVmod-ZV = 0). (ii) Every fermionic hole 
(open circle) inside represents a spin-up spinon (square surrounding open circle) and 
every fermionic particle (full circle) outside represents a spin-down spinon (square 
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Figure Al. Fermion configurations of all eigenstates for N = 4 of the XX model 
(J z = 0). Fermionic particles (holes) are denoted by full (open) circles. Spinons 
with spin up (down) are denoted by (black or gray) squares around open (full) 
circles. Solitons with spin up (down) are denoted by black squares around open 
(full) circles. The fermion momenta rhi (in units of tt/N) can be read off the 
diagram. The spinon orbital momenta rrii (also in units of tt/N) and the spinon 
spins &i are given explicitly and can be inferred from the fermion configuration 
as explained in the text. Also given are the wave number k (in units of 2n/N) 
and the energy E (in units of J±) of each eigenstate. 
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surrounding full circle), (iii) Any number of adjacent spinons in the representation of 
Fig. IA1I are in the same orbital. Two spin-up (spin-down) spinons that are separated 
by £ consecutive fermionic particles (holes) have quantum numbers separated by 11. 
(iv) The spinon orbital momenta are sorted in increasing order from the right-hand 
prong of the gray fork toward the left across the inside domain and toward the right 
with wrap-around through the outside domain. 

For the determination of the soliton content of any XX eigenstate we must also 
consider the black fork in Fig. IA11 which is shifted two units to the left. If there is 
any spin- up/spin-down spinon pair as identified by the rules pertaining to the gray 
fork that does not also satisfy the same rules with respect to the black fork, then it 
is omitted from the list of solitons. All such spinons are identified by gray squares 
m Fig. EU The spinon vacuum is just one of two soliton vacua, the other soliton 
vacuum being a two-spinon state. The two soliton vacua have wave numbers differing 
by 7r. In the Ising limit of the XXZ model they correspond to the symmetric and 
antisymmetric combinations of the two product Neel states. 
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